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Abstract 

Recent HI data on triple differential dijet cross sections in e^p interactions in the region of 
low photon virtualities are shown to be in reasonable agreement with the predictions of the 
NLO QCD calculations obtained using the program NLOJET++. The implications of this 
observation for the phenomenological relevance of the concept of resolved virtual photon 
are discussed. 



1 Introduction 



Recent HI data on triple differential dijet cross sections in e ± p interactions in the region of 
low photon virtualities [ 1 1 have revealed a clear excess of the data over NLO QCD predictions 
obtained using the DISENT program [2|. This excess comes predominately from the region of 
low photon virtualities Q 2 (Q 2 > 2 GeV 2 ) and small x 7 (x 7 < 0.75), where x 7 denotes fraction 
of the four- momentum of the photon carried by the parton involved in the hard collision 1 . We 
recall that DISENT dispenses with the concept of resolved photon and evaluates, up the order 
aa 2 , the direct photon contribution only. In the region Q 2 > 2 GeV 2 this is quite a legitimate 
procedure. However, it has been argued EJ that even for such moderate values of Q 2 the 
concept of resolved virtual photon is useful phenonomenologically as a way of resuming part 
of higher order direct photon contributions. These higher order QCD corrections are impor- 
tant particularly in the region s 7 < 0.75, to which only the lowest order tree level diagrams 
contribute in DISENT. 

In this region DISENT calculations are described by diagrams exemplified by those in Fig. 
[T^,c. In the kinematic region t\ <C r 2 ~ E\ (where Ti, r 2 denote the virtualities of the ladder 
partons in Fig. [Tk,c) the contribution to the dijet cross section coming from the diagram in Fig. 
[T^ can be approximated by the sum of the convolutions which describe the contributions of 
transverse (k = T) and and longitudinal (k = L) polarizations of the virtual photon. 
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describe the fluxes of transverse and longitudinal photons in the incoming electron, the QED 
contributions to the quark and gluon distribution functions of the transverse and longitudinal 
virtual photons are given as 
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and <j LO (qq — > GG) stands for the lowest order contribution to the cross section of the process 
qq —> GG. Similar expressions can be written for other partonic subprocesses. In Eqs. <El|6l), 
ei denotes the electric charge of the quark qi and M stands for the factorization scale, in ([TJ 
identified with jet transverse energy. The full expressions for the distribution functions <|4l|6l), 
containing the exact Q 2 dependence with the correct threshold behaviour for Q 2 /m 2 — > 0, can 
be found in H. The contributions of the diagrams in Fig. [Tk,c contain also singularities coming 
from the region of small r 2 and those of Fig. Eb,d from the region of small r 3 , but these are 
understood to be absorbed in the parton distribution functions of the proton. 

'We use the notation and terminology as described in detail in Q. 
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Figure 1: Examples of the LO (a) and (c) and NLO (b) and (d) Feynman diagrams included 
in NLOJET++ and contributing to the three jet cross section in the region x 7 < 0.75. In the 
framework of the concept of resolved virtual photon the integrals over the region of small n 
give rise to quark (a-c) and gluon (d) distribution function of the photon. 

In ID it has also been shown that adding in HERWIG the LO resolved photon contribution 
to the LO direct one helps bring the QCD predictions closer to the data. At the NLO such 
a procedure can be performed only at the parton level and currently only within the JETVIP 
program [ 5 1 . The NLO QCD corrections come loop corrections to LO diagrams like those of 
Fig. [Tk,c and from tree Feynman diagrams exemplified by those of Fig. Ch,d, which contain 
terms proportional to the "large collinear logs" of the form In 2 [E\ /r 2 ). Within the framework of 
resolved virtual photon contribution, these terms can be interpreted as lowest QCD corrections 
to the purely QED part of quark and gluon distribution functions of the photon (EE). In Fig. 
El they correspond to the second (first) term in the definition of the so called pointlike part of 
quark (gluon) distribution function (D^Jl ) of the photon Q. 

Contrary to DISENT, which uses the subtraction method, JETVIP employs the phase space 
slicing method to regularize mass singularities. Unfortunately, the NLO resolved photon cal- 
culations obtained with JETVIP turned out [6 1 to be rather sensitive to choice the associated 
slicing parameter y c . Consequently, although adding resolved photon contribution worked in 
the right direction, no quantitative conclusion could be drawn (Q. 

The situation has recently changed due to the appearance of the NLOJET++ program 0, 
which allows the user to calculate beside single and dijet cross sections also the triple jet cross 
sections to the NLO. This program uses the same subtraction method as DISENT and similarly 
as the latter dispenses with the concept of resolved virtual photon. As far as single and double 
jet cross sections are concerned it is thus in principle identical to DISENT. It, however, provides 
also the option - which we shall call 3 -jet mode - of calculating NLO QCD correction to three 
well separated jets. This mode can be selected by cutting off the region of x-y close to 1, where 
also the dijet final states contribute. Respecting the binning of the data in HI we set x 7 < 0.75. 
Comparing the results obtained with NLOJET++ and JETVIP we could address two questions: 

• how important are terms of the order aa^ included in exact NLO calculations, 

• how important is the all order resummation of the collinear logarithms into the QCD- 
improved PDF of the photon used in the resolved photon contribution of JETVIP? 

Unfortunately, as mentioned above the NLO resolved photon contribution are rather sensitive to 
the variation of the slicing parameter y c . Moreover, there is also a sizable discrepancy between 
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Figure 2: Graphs contributing to the definition of the pointlike part of quark and gluon distribu- 
tion functions of the photon. 

the results of the NLO direct photon contribution of JETVIP on one side and those of DISENT 
or NLOJET++ in 2-jet mode on the other, with JETVIP results lying in the region x 7 < 0.75 
significantly below those of DISENT or NLOJET++. This discrepancy, clearly visible in Fig. 
is not understood. 



2 Data Selection 

The data employed in this article have been published in JT), where a detailed description of 
event selection and experimental methods used in their extraction can be found. Only the main 
facts are therefore recalled. The data, taken at HERA in the years 1999 and 2000, when elec- 
trons with an energy of 27.55 GeV collided with protons with an energy of 920 GeV, correspond 
to an integrated luminosity of 57 pb _1 . The kinematic region of the present analysis is defined 
by cuts on the photon virtuality Q 2 and its inelasticity y 

2 < Q 2 < 80 GeV 2 , 0.1 < y < 0.85 (7) 

as well as by cuts on the hadronic final state, which has to contain at least two jets found 
using the longitudinally invariant k t jet algorithm 0. The jet transverse energies, E^ , and 
pseudorapidities, rj* , are calculated relative to the 7*p collision axis in the YP center-of-mass 
frame. The jets are ordered according to their transverse energy, with jet 1 being the highest 
jet. The two jets with the highest transverse energies (leading jets) are required to satisfy 

E* T1 >7 GeV, E* T2 > 5 GeV, (8) 



-2.5 < tjI < 0, - 2.5 < 7/ 2 * < 0. (9) 
In total 105 658 events satisfied these selection criteria. For dijet events the variable 

E (E* ~ Pi,,) 

E ( E -Pi) 

hadrons 

was used as a hadron level estimate of x 7 . The sum in the numerator runs over the two leading 
jets and the sum in the denominator includes the full hadronic final state. 

The data were corrected for initial and final state QED radiation effects, trigger inefficien- 
cies, limited detector acceptance and resolution and a photoproduction background. The sys- 
tematic errors of different origin are added in quadrature. The dominant source of the systematic 
error arises from the uncertainty of the energy calibration of the HI calorimeters and the model 
dependence of acceptance corrections Q. 
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3 Results 



In Fig. |3]the NLOJET++ results corresponding to both 2-jet and 3-jet modes and evaluated 
in the kinematic region (J?© are compared with the HI data [1|. For comparison, the results 
obtained with DISENT [ 1 1 are plotted as well. 

There is a large difference between the NLOJET++ results obtained in 2-jet and 3-jet modes, 
the latter lying systematically above the former. This difference, which is most pronounced for 
small x 7 and low Q 2 , indicates the importance of the NLO QCD corrections in this region. The 
NLOJET++ results in the 3-jet mode come significantly closer to the data than those of the 
2-jet mode, but they still undershoot it. The remaining gap between the data and NLOJET++ 
calculations in the 3-jet mode is again most pronounced for small x 7 and low Q 2 . In view of 
large NLO QCD corrections in this region this remaining discrepancy is not surprising and and 
allows room for still higher order QCD corrections. 

At the moment the only, though approximate, way of estimating these higher order terms 
exploits the concept of resolved virtual photon and the only code that offers this option at the 
NLO is JETVIP. Despite the problems mentioned above, we plot in Fig. |3] the full JETVIP 
results Q, i-e- the sum of direct (with photon splitting term subtracted) and resolved photon 
contributions. Ideally, one would expect the full JETVIP results to be above those of NLO- 
JET++ in 3-jet mode. Fig. |3]shows that they are in fact close to each other. This may be caused 
by various effects, but we would like to point out the following one. As already mentioned 
above the direct photon results of JETVIP are systematically below those of NLOJET++ in the 
two jet mode. Assuming the latter are correct but the resolved photon contribution is estimated 
correctly by JETVIP, the "correct" full JETVIP should lie somewhat above the NLOJET++ 3- 
jet mode predictions and close to the data. At the moment, this is just a pure speculation, but 
we hope to check it, once a new version of JETVIP becomes available [ 9 1 . 

In summary, we have shown that the NLOJET++ calculations of the dijet cross sections in 
the 3-jet mode are significantly closer to the HI data 1 1 1 that those of DISENT. This demon- 
strates the importance of the NLO QCD corrections in processes involving virtual photons. The 
remaining gap between the data and current QCD calculations in the region of low Q 2 and small 
x y may be closed using the concept of the resolved virtual photon contribution. 
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Figure 3: Triple differential dijet cross section d 3 o"2jet/dQ 2 d£'y dx J 7 ets with asymmetric cuts 
(see text). The inner error bars on the data points show the statistical error, the outer ones the 
quadratic sum of systematic and statistical errors. The data are compared to NLO direct photon 
calculations using DISENT (full line) and JETVIP (dotted line), the sum of NLO direct and 
resolved photon contributions of JETVIP (dashed line) and the NLOJET++ predictions in 2-jet 
(triangles) and 3-jet (squares) mode. All calculations are corrected for hadronisation effects. 
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